IMMUNOLOGY2013” Travel Grants 


tom R&D Systems . ., ++ 


Comparative Analysis of Bat Genomes Provides Insight into the 
Evolution of Flight and Immunity 

Guojie Zhang et al. 

Science 339, 456 (2013); 

DOI: 10.1126/science.1230835 


This copy is for your personal, non-commercial use only. 


If you wish to distribute this article to others, you can order high-quality copies for your 
colleagues, clients, or customers by clicking here. 


Permission to republish or repurpose articles or portions of articles can be obtained by 
following the guidelines here. 


The following resources related to this article are available online at 
www.sciencemag.org (this information is current as of January 24, 2013 ): 


Updated information and services, including high-resolution figures, can be found in the online 
version of this article at: 
http:/www.sciencemag.org/content/339/61 18/456. full.html 


Supporting Online Material can be found at: 
http://www.sciencemag.org/content/suppl/201 2/1 2/1 9/science. 1230835.DC1.html 


This article cites 51 articles, 28 of which can be accessed free: 
http://www.sciencemag.org/content/339/61 18/456. full.html#ref-list-1 


This article appears in the following subject collections: 
Genetics 
http://www.sciencemag.org/cgi/collection/genetics 


Science (print ISSN 0036-8075; online ISSN 1095-9203) is published weekly, except the last week in December, by the 
American Association for the Advancement of Science, 1200 New York Avenue NW, Washington, DC 20005. Copyright 


2013 by the American Association for the Advancement of Science; all rights reserved. The title Science is a 
registered trademark of AAAS. 


Downloaded from www.sciencemag.org on January 24, 2013 


REPORTS 


456 


giving the axonal cytoskeleton a long-range 
order. Despite the molecular composition dif- 
ferences between the axon initial segments and 
distal axons [for example, ankyrin-G and BIV- 
spectrin are confined in the axon initial segment 
by an exclusion effect of the distal axon proteins 
ankyrin-B and fll-spectrin (3/)], the cytoskel- 
etal organization is similar between the initial 
and distal segments of the axons, both adopting 
a quasi-1D, periodic structure. Interestingly, we 
found this periodic cytoskeleton structure to be 
present only in axons, not in dendrites, which 
instead primarily contained long actin filaments 
running along the dendritic axis. Although the 
microscopic interactions between the molecu- 
lar components of the axon cytoskeleton are 
probably similar to those between the erythro- 
cyte analogs (9, 10), the overall structure of this 
quasi-1D, periodic cytoskeleton in axons is 
distinct from the 2D, pentagonal or hexagonal 
structure observed for the erythrocyte membrane 
cytoskeleton (//, 12). In Drosophila motoneuron 
axons near the neuromuscular junctions, spectrin 
and ankyrin appear to organize into an erythrocyte- 
like, pentagonal or hexagonal lattice structure (/6), 
which is distinct from the quasi-1D, periodic, lad- 
derlike structure that we observed in the axons of 
vertebrate brains. Whether the difference is due 
to invertebrate versus vertebrate animals or pe- 
ripheral versus central nervous systems is a topic 
for future investigations. 

The periodic, actin-spectrin—based cytoskel- 
eton observed here may not be involved in myosin- 
dependent axonal transport. If the analogy to the 
erythrocyte membrane cytoskeleton holds, the 
capped short actin filaments in the ringlike actin 
structures in axons are probably bound by tropo- 
myosin (9, 10), which could potentially prevent 
the binding of myosins. Myosin-dependent ax- 
onal transport could, however, be mediated by 
the long actin filaments that run along the axon 
shaft. The quasi-1D, periodic, actin-spectrin cyto- 
skeleton may instead provide elastic and stable 
mechanical support for the axon membrane, giv- 
en the flexibility of spectrin. Elastic and stable 
support is particularly important for axons, be- 
cause they can be extremely long and thin and 
have to withstand mechanical strains as ani- 
mals move (37). Indeed, the loss of B-spectrin 
in Caenorhabditis elegans leads to spontaneous 
breaking of axons, which is caused by mechan- 
ical strains generated by animal movement and 
can be prevented by paralyzing the animal (37). 
The highly periodical submembrane cytoskeleton 
can also influence the molecular organization of 
the plasma membrane by organizing important 
membrane proteins along the axon. We found 
that sodium channels were distributed periodically 
along the axon initial segment in a coordinated 
manner with the underlying actin-spectrin cyto- 
skeleton. An axonal plasma membrane with pe- 
riodically varying biochemical and mechanical 
properties may not only influence how an action 
potential is generated and propagated, but might 
also affect how axons interact with other cells. 
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Bats are the only mammals capable of sustained flight and are notorious reservoir hosts 

for some of the world’s most highly pathogenic viruses, including Nipah, Hendra, Ebola, and 
severe acute respiratory syndrome (SARS). To identify genetic changes associated with the 
development of bat-specific traits, we performed whole-genome sequencing and comparative 
analyses of two distantly related species, fruit bat Pteropus alecto and insectivorous bat 

Myotis davidii. We discovered an unexpected concentration of positively selected genes in the 
DNA damage checkpoint and nuclear factor xB pathways that may be related to the origin of flight, 
as well as expansion and contraction of important gene families. Comparison of bat 

genomes with other mammalian species has provided new insights into bat biology and evolution. 


ats belong to the order Chiroptera with- 
B in the mammalian clade Laurasiatheria 

(1). Although consensus has not been 
reached on the exact arrangement of groups with- 
in Laurasiatheria, a recent study placed Chiroptera 
as a sister taxon to Cetartiodactyla (whales + even- 
toed ungulates such as cattle, sheep, and pigs) 
(2). The Black flying fox (Pteropus alecto) and 
David’s Myotis (Myotis davidii) represent the 
Yinpterochiroptera and Yangochiroptera subor- 


ders, respectively, and display a diverse range of 
phenotypes (Fig. 1). Captive colonies, immor- 
talized cell lines, and bat-specific reagents have 
been developed for these two species; however, 
genomic data are currently unavailable. 

The most conspicuous feature of bats, distin- 
guishing them from all other mammalian species, 
is the capacity for sustained flight. Positive selec- 
tion in the oxidative phosphorylation (OXPHOS) 
pathway suggests that increased metabolic capac- 
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ity played a key role in its evolution (3), yet the 
by-products of oxidative metabolism [such as re- 
active oxygen species (ROS)] can produce harm- 
ful side effects including DNA damage (4). We 
hypothesize that genetic changes during the evolu- 
tion of flight in bats likely included adaptations 
to limit collateral damage caused by by-products 
of elevated metabolic rate. Another phenome- 
non that has sparked intense interest in recent 
years is the discovery that bats maintain and dis- 
seminate numerous deadly viruses (5). In this con- 
text, we further hypothesize that the long-term 
coexistence of bats and viruses must have im- 
posed strong selective pressures on the bat ge- 
nome, and the genes most likely to reflect this 
are those directly related to the first line of anti- 
viral defense—the innate immune system. 

We performed high-throughput whole-genome 
sequencing of individual wild-caught specimens 
of P. alecto and M. davidii using the Illumina 
HiSeq platform (6). More than 100 x coverage 
high-quality reads were obtained for P. alecto 
and M. davidii, which resulted in high-quality 
assemblies (tables S1 to S3 and fig. S1). The two 
bat genomes, at ~2 Gb, were smaller in size than 
other mammals (7) (fig. S2), whereas the number 
of genes we identified was similar to those of other 
mammals (21,392 and 21,705 in P. alecto and 
M. davidii, respectively) (fig. S3). Both species 
displayed a high degree of heterozygosity at the 
whole-genome level (0.45% and 0.28% in P alecto 
and M. davidii, respectively) (tables S4 and S5), 
whereas repetitive content accounted for slight- 
ly less than one-third of each genome (tables S6 
and S7). We identified a novel endogenous viral 
element derived from Saimiriine herpesvirus 2 
that has expanded to 126 copies in P. alecto 
(table S8 and fig. S4). Gene family expansion and 
contraction analysis (tables S9 to S12) revealed 
significant expansion (P < 0.05) of 71 gene fam- 
ilies in M. davidii compared with only 13 in P alecto, 
which may be related to a recent wave of DNA 
transposon activity (8). 

We screened all nuclear-encoded bat genes to 
identify those for which a single orthologous copy 
was unambiguously present in both bat species 
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Fig. 1. Comparison of bat biological traits. P. alecto and M. davidii represent two distinct Chiropteran 
suborders and demonstrate diverse evolutionary adaptations. PNG, Papua New Guinea. 
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Fig. 2. Phylogenomic analysis. Maximum -likelihood phylogenomic analysis of 2492 genes from M. davidii, 
P. alecto, and eight mammalian species. Divergence time estimates in blue, gene family expansion 
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as well as in human, rhesus macaque, mouse, 
rat, dog, cat, cattle, and horse. From this, 2492 
genes were used to perform maximum-likelihood 
and Bayesian phylogenomic analysis (Fig. 2 and 
figs. S5 to S7). All phylogenetically informative 
signals, including concatenated nucleotides and 
amino acids, vigorously supported bats as a mem- 
ber of Pegasoferae (Chiroptera + Perissodactyla + 
Carnivora) (9), with the bat lineage diverging from 
the Equus (horse) lineage ~88 million years ago, 
buttressed by findings at the transcript level (/0). 
However, phylogenetic reconstruction with mito- 
chondrial DNA sequences resulted in bats oc- 
cupying an outlying position in Laurasiatheria 
(fig. S8). The incongruence between nuclear and 
mitochondrial trees likely reflects rapid evolution 
of the mitochondrial genome of the bat ancestor 
during the evolution of flight (3). 

To identify mechanisms that facilitated the 
origin of flight in bats, we surveyed genes in- 
volved in detection and repair of genetic damage. 
A high proportion of genes in the DNA damage 
checkpoint-DNA repair pathway were found to 
be under positive selection in the bat ancestor, 
including ATM, the catalytic subunit of DNA- 
dependent protein kinase (DNA-PKc), RAD5O, 
KUS80, and MDM2 (Fig. 3A and Table 1). We 
propose that these changes may be directly re- 
lated to minimizing and/or repairing the negative 
effects of ROS generated as a consequence of 
flight. Additionally in this pathway, TP53 (p53) 
and BRCA2 were shown to be under positive 
selection in M. davidii, whereas LIG4 was under 
positive selection in P. alecto (Table 1). Bat- 
specific mutations in a nuclear localization signal 
in p53 and a nuclear export signal in MDM2 (Fig. 
3B and fig. S9) may affect subcellular localiza- 
tion and function in both species (//, 12). Other 
candidate flight-related genes under positive se- 
lection in the bat ancestor included COL3A/, in- 
volved in skin elasticity, and CACNA2D1, which 
has a role in muscle contraction (table $13). 

We next examined genes of the innate im- 
mune system (Table 1). Positively selected genes 
in the bat ancestor included c-REL, a member 
of the nuclear factor kB (NF-«B) family of tran- 
scription factors, which also contained amino acid 
changes potentially affecting inhibitor of NF-KB 
(IkB) binding (fig. S10). In addition to diverse 
roles in innate and adaptive immunity (13), c-REL 
plays a role in the DNA damage response by 
activating ATM (14) and CLSPN (15), whereas 
ATM is also an upstream regulator of NF-«B 
(16). The DNA damage response plays an im- 
portant role in host defense and is a known target 
for virus interaction (/7), which raises the pos- 
sibility that changes in DNA damage response 
mechanisms during selection for flight could 
have influenced the bat immune system. 

It is intriguing that both P alecto and M. davidii 
have lost the entire locus containing the PYHIN 
gene family, including AJM2 and JFI/6, both of 
which are involved in sensing microbial DNA 
and the formation of inflammasomes (fig. S11). 
The association between PYHIN genes and cell 
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Fig. 3. Accelerated evolution in the DNA damage checkpoint in bats. (A) Positive selection in the DNA 
damage checkpoint—DNA repair pathway. Genes under positive selection in the bat ancestor are high- 
lighted in orange. Genes under positive selection in M. davidii only (p53, BRCA2) or P. alecto only (LIG4) 
are highlighted in blue. IFN, interferon; IL, interleukin. (B) Mutations unique to bats were detected in the 
functionally relevant regions of the p53 nuclear localization signal (NLS) and MDM2 nuclear export 


signal (NES) (black shading). 


cycle regulation in other species (/8) hints that 
loss of the PYHIN family in bats may be con- 
nected to changes in the DNA damage pathway, 
because at least one PYHIN gene is present in 
all other major groups of eutherian mammals 
(19). NLRP3, triggered by both viral infection and 
ROS in other mammals (20), plays an analogous 
role to AJM2 in inflammasome assembly and was 
also under positive selection in the bat ancestor 
(Table 1). 

Natural killer (NK) cells provide a first line of 
defense against viruses and tumors and include 
two families of NK cell receptors; killer-cell 
immunoglobulin like receptors (KIRs), encoded 
by genes in the leukocyte receptor complex 


(LRC), and killer cell lectin-like receptors (KLRs, 
also known as Ly49 receptors), encoded within 
the natural killer gene complex (NKC). KLRs 
and KIRs were entirely absent in P. alecto and 
reduced to a single Ly49 pseudogene in M. davidii 
(table S14). KIR-like receptors identified in other 
species (2/) were also absent from both P. alecto 
and M. davidii genomes, which was supported 
by transcript analysis in P alecto (10). This like- 
ly indicates that bat NK cells use a novel class of 
receptors to recognize classical major histocom- 
patibility complex class I molecules. Furthermore, 
additional LRC members of the immunoglobulin 
superfamily [including sialic acid-binding im- 
munoglobulin-like lectins (SIGLECs), leukocyte 


25 JANUARY 2013 VOL 339 SCIENCE www.sciencemag.org 


Downloaded from www.sciencemag.org on January 24, 2013 


REPORTS 


cattle, and horse. œ0 is the average ratio in all branches, œ1 is the 
average ratio in nonbat branches, and 2 is the ratio in the bat branch. A 
low P value indicates that the w2 model fits the data better than the 


Table 1. DNA damage checkpoint and innate immune genes under pos- 
itive selection in the bat lineages. The rate ratio œ of dN/dS was calculated 
using multiprotein alignments of P. alecto and M. davidii sequences with 


orthologous sequences from human, rhesus macaque, mouse, rat, dog, 1 model. 
; o0 o1 a2 
Lineage Symbol Gene Gverage) lothar) (target) P value 
Ancestor TLR7 Toll-like receptor 7 0.2821 0.2670 2.7778 3.54E-07 
ATM Ataxia telangiectasia mutated 0.20096 0.19595 0.7163 1.34E-05 
MDM2 Mdm2 p53 binding protein homolog (mouse) 0.13358 0.12615 0.81085 4.05E-04 
NLRP3 NLR family, pyrin domain—containing 3 0.1788 0.1714 1.1884 1.93E-04 
MAP3K7 Mitogen-activated protein kinase kinase kinase 7 0.0216 0.0194 0.4786 8.93E-03 
RAD50 RAD50 homolog 0.09657 0.09343 0.28882 7.95E-03 
PRKDC Protein kinase, DNA-activated, catalytic polypeptide 0.23036 0.22768 0.45155 6.80E-03 
KU80 X-ray repair complementing defective repair in Chinese hamster cells 5 0.31145 0.30436 0.91747 3.75E-02 
c-REL v-rel reticuloendotheliosis viral oncogene homolog (avian) 0.2495 0.2403 1.5717 1.11E-02 
P. alecto TBK1 TANK-binding kinase 1 0.0643 0.0522 0.2930 1.29E-09 
LIG4 Ligase IV, DNA, ATP-dependent 0.12033 0.11376 0.24797 8.91E-04 
IL18 Interleukin 18 (interferon-y—inducing factor) 0.5298 0.4532 1.7647 2.66E-04 
IFNG Interferon-y 0.5010 0.4527 1.3282 4.89E-03 
ISG15 ISG15 ubiquitin-like modifier 0.2069 0.1909 0.4387 2.63E-02 
DDX58 DEAD (Asp-Glu-Ala-Asp) box polypeptide 58 0.3040 0.2923 0.4661 1.23E-02 
M. davidii IFNAR1 Interferon (o, B, and œ) receptor 1 0.4954 0.4723 31.0924 7.00E-03 
TP53 Tumor protein p53 0.25623 0.23933 0.48123 7.00E-03 
BRCA2 Breast cancer 2, early onset 0.49002 0.47732 0.64213 1.31E-03 
IRAK4 Interleukin-1 receptor—associated kinase 4 0.1670 0.1583 0.3531 1.96E-02 


immunoglobulin-like receptors (LILRs), carcino- 


lection in M. davidii and three other hibernating 


3. Y. Y. Shen et al., Proc. Natl. Acad. Sci. U.S.A. 107, 


embryonic antigen-related cell adhesion mol- species relative to nonhibernators (table S16). 8666 (2010). : . 
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dergone considerable gene duplication in M. davidii AP-1 transcription factor), had significantly high- 649 (2011). 


and other mammals yet have almost complete- 
ly failed to expand in P. alecto (fig. S12). As the 
genes encoded within the LRC bind a variety of 
ligands and play multiple roles in immune regu- 
lation, these observations have diverse implica- 
tions for differences in immune function between 


er ratio of nonsynonymous to synonymous sub- 
stitutions (dN/dS) in the echolocating M. davidii 
branch relative to non-echolocating branches 
(table S17). Of note, the third exon in M. davidii 
FOXP2 had even greater variation from the mam- 
malian consensus than two previously identified 


6. Materials and methods are available as supplementary 
materials on Science Online. 
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RNASE4 genes (re ative to the mammalian con- role in the evo ution o Myotis bats and may have 17. A. S. Turnell, R. J. Grand, J. Gen. Virol. 93, 2076 
sensus) that suggest diversification of substrate helped contribute to their speciation. Concen- (2012). 
specificity (fig. S13). With a proven roleinhost tration of positively selected genes in the DNA 18. S$. A. Schattgen, K. A. Fitzgerald, Immunol. Rev. 243, 
: : 2 : i aadi 109 (2011). 
defense against RNA viruses (22), RNASE4 ex: damage checkpoint pathway in bats may indi 19. 1. A. Cridland et al., BMC Evol. Biol 12, 140 (2012), 
pansion in M. davidii may have implications for cate an important step in the evolution of flight, 20, J. Tschopp, K. Schroder, Nat. Rev. Immunol. 10, 210 
virus resistance but may also reflect the insec- whereas evidence of change in components shared (2010). 
tivorous diet of M. davidii, in contrast with that by the DNA damage pathway and the innate im- 21- i G. T S. Beck, Curr. Opin. Immunol. 19, 553 
: : j : : : eee 2007). 
of P. alecto, which consumes predominantly fruit, mune system raises the interesting possibility that 55 F. Cocchi et al., Proc. Natl. Acad. Sci. U.S.A. 109, 5411 
flowers, and nectar. flight-induced adaptations have had inadvertent (2012). 
M. davidii also differs from P. alecto in aspects effects on bat immune function and possibly also 23. G. Li, J. H. Wang, S. J. Rossiter, G. Jones, S. Y. Zhang, 
including hibernation and echolocation (Fig. 1). life expectancy (24). The data generated by this PLoS ONE 2, e900 (2007). l 
24. A. K. Brunet-Rossinni, S. N. Austad, Biogerontology 5, 


Bile salt-stimulated lipase (BSSL), capable of 
hydrolyzing triglycerides into monoglycerides 
and subsequently releasing digestible free fatty 
acids, has been specifically expanded in M. davidii 
compared with P alecto and other mammals (fig. 
S14). In addition, we observed six candidate genes 
related to hibernation, which showed positive se- 


study will help to address major gaps in our un- 
derstanding of bat biology and to provide new 
directions for future research. 


References and Notes 

1. E. C. Teeling et al., Science 307, 580 (2005). 

2. M. F. Nery, D. J. González, F. G. Hoffmann, J. C. Opazo, 
Mol. Phylogenet. Evol. 64, 685 (2012). 


211 (2004). 


Acknowledgments: We thank H. Field, C. Smith, and M. Yu 
for helping source genomic DNA; K. Itahana and J. J. Boomsma 
for constructive discussion; and M. Cowled for graphics 
assistance. We acknowledge financial support from the China 
National Genebank at Shenzhen, CSIRO (Office of the Chief 
Executive Science Leaders Award, Julius Award), The Australian 


www.sciencemag.org SCIENCE VOL 339 25 JANUARY 2013 


459 


Downloaded from www.sciencemag.org on January 24, 2013 


REPORTS 


460 


Research Council (FT110100234), State Key Program for 
Basic Research (2011CB504701), National Natural Science 
Foundation of China (81290341), and the Defense Threat 
Reduction Agency of the USA. The views expressed in this 
article are those of the authors and do not necessarily reflect 
the official policy or position of the Department of the Navy, 
Department of Defense, or the U.S. government. K.A.B.-L. 
and K.G.F. are contractors for the U.S. government. This work 
was prepared as part of their official duties. Title 17 U.S.C. 
§105 provides that “Copyright protection under this title is 
not available for any work of the United States Government.” 
Title 17 U.S.C. §101 defines a U.S. government work as a work 
prepared by a military service member or employee of the 
U.S. government as part of that person’s official duties. 

P. alecto and M. davidii genomes have been deposited at 
DNA Data Bank of Japan/European Molecular Biology 
Laboratory/GenBank under the accession nos. ALWS01000000 


and ALWT01000000. Short-read data have been deposited 
into the Short Read Archive under accession nos. SRA056924 
and SRA056925. Raw transcriptome data have been deposited 
in Gene Expression Omnibus as GSE39933. Tree files and 
alignments have been submitted to TreeBASE under Study 
Accession URL: http://purl.org/phylo/treebase/phylows/study/ 
1B2:S13654. We also thank the editor and two anonymous 
reviewers for their helpful comments and suggestions. 
Author contributions: J.W., L.-F.W., G.Z., C.C.B., and K.A.B.-L. 
conceived the study. M.T., M.L.B., G.A.M., G.C., L.W., and 
Z.S. prepared the samples. G.Z., Z.H., X.F., Z.X., W.Z., Y. Zhu, 
X.J., LY., J.X. Y.F., Y.C., X.S., Y. Zhang, K.G.F., K.A.B.-L., 
and J.W. performed genome sequencing, assembly, and 
annotation. G.Z. and J.W. supervised genome sequencing, 
assembly, and annotation. G.Z., C.C., Z.H., X.F., J.W.W., Z.X., 
J.N., W.Z., P.Z., Y. Zhu, M.T., and M.L.B. performed genome 
analyses. G.Z., Z.H., C.C., and J.W.W. carried out genetic 


analyses. G.Z., C.C., Z.H., X.F., P.Z., J.N., M.T., J.W.W., M.L.B., 
and L.-F.W. discussed the data. All authors contributed to data 
interpretation. C.C. and J.W.W. wrote the paper with 
significant contributions from G.Z., Z.H., P.Z., J.N., M.T., 
M.L.B., and L.-F.W. and input from all authors. The authors 
declare no competing financial interests. Requests for materials 
should be addressed to the authors for correspondence. 


Supplementary Materials 
www.sciencemag.org/cgi/content/full/science.1230835/DC1 
Materials and Methods 

Figs. $1 to $15 

Tables S1 to $17 

References (25-52) 


28 September 2012; accepted 10 December 2012 
Published online 20 December 2012; 
10.1126/science.1230835 


Tunable Signal Processing Through 
Modular Control of Transcription 
Factor Translocation 


Nan Hao,”? Bogdan A. Budnik,* Jeremy Gunawardena,’ Erin K. O’Shea*2* 


Signaling pathways can induce different dynamics of transcription factor (TF) activation. We explored 
how TFs process signaling inputs to generate diverse dynamic responses. The budding yeast general 
stress—responsive TF Msn2 acted as a tunable signal processor that could track, filter, or integrate 
signals in an input-dependent manner. This tunable signal processing appears to originate from dual 
regulation of both nuclear import and export by phosphorylation, as mutants with one form of 
regulation sustained only one signal-processing function. Versatile signal processing by Msn2 is 
crucial for generating distinct dynamic responses to different natural stresses. Our findings reveal how 
complex signal-processing functions are integrated into a single molecule and provide a guide for 


input amplitude was chosen on the basis of the 
steady-state amount of Msn2 nuclear localization 
in response to sustained inputs: high-amplitude 
input (3 uM 1-NM-PP1) led to maximal nuclear 
localization of Msn2, whereas low-amplitude 
input (0.2 uM 1-NM-PP1) induced an inter- 
mediate amount of nuclear localization (Fig. 1B, 
black circles). The pulse duration of oscillatory 
input was selected on the basis of duration of 
pulsatile Msn2 nuclear bursts in the physiolog- 
ical response to glucose limitation (6). With high- 
amplitude oscillatory input, each input pulse 
induced a large amount of nuclear localization 
(Fig. 1C, left). In contrast, oscillatory input with 
low amplitude barely elicited any localization 
responses, although sustained input with the same 
amplitude led to a half-maximal amount of nu- 


the design of TFs with “programmable” signal-processing functions. 


any transcription factors (TFs) display 
Mi= activation dynamics in response 
to various external stimuli (7—4). To in- 
vestigate how TFs process upstream signals, we 
studied the Saccharomyces cerevisiae general 
stress-responsive TF Msn2 (5). In the absence 
of stress, Msn2 is phosphorylated by protein 
kinase A (PKA) and localized to the cytoplasm; 
in response to stress, Msn2 is dephosphorylated 
and translocates to the nucleus, where it induces 
gene expression (5). 
Natural stresses elicit highly variable dynam- 
ics of Msn2 nuclear translocation (Fig. 1A) (6, 7), 
which are thought to result from oscillatory sig- 
naling inputs (presumably PKA activity) (8). To 
study how Msn2 processes oscillatory PKA in- 
puts, we used an engineered yeast strain (6) carry- 
ing mutations in all three PKA isoforms that 
enable selective inhibition of PKA activity by a 
cell-permeable inhibitor, 1-NM-PP1 (9). We used 
this synthetic system and a microfluidics plat- 
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clear localization (Fig. 1C, right). Therefore, Msn2 
filters temporal fluctuations of the input in an 
amplitude-dependent manner such that it tracks 
high-amplitude inputs, but responds in a limited 
manner to low-amplitude signals. 


form (/0) mounted on a microscope to produce 
oscillatory inputs of PKA inhibition and moni- 
tored translocation of Msn2 to the nucleus. The 
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Fig. 1. Tunable signal-processing behaviors of Msn2. (A) Illustration of 
the distinct single-cell dynamic responses of Msn2 to various stresses. (B) 
Steady-state abundance of Msn2 in the nucleus in response to various 
concentrations of 1-NM-PP1. In response to each concentration of 1-NM-PP1, Msn2 exhibited uniform 
and stable nuclear localization in single cells and did not exhibit stochastic fluctuations as observed in 
natural stress responses. Open circles: responses to different concentrations of 1-NM-PP1; closed circles: 
responses to 3 uM and 0.2 uM 1-NM-PP1, which are used as high- and low-amplitude inputs, respectively, 
for the following analyses. AU, arbitrary unit. (C) Averaged single-cell time traces of Msn2 nuclear trans- 
location (bottom: n = 50 cells; error bar: single-cell variances) in response to oscillatory inputs with high 
and low amplitudes (top). (Left) High-amplitude input produced by 3 uM 1-NM-PP1; (right) low-amplitude 
input produced by 0.2 uM 1-NM-PP1. Pulse duration of 3 min; pulse interval of 2 min. To emphasize the 
fact that 3 uM 1-NM-PP1 elicits a steady-state response that is about twice the response elicited by 0.2 uM 
1-NM-PP1, the top y axes are not presented on a linear scale. 
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